Three-dimensional angular momentum projection in relativistic mean-field theory 
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Based on a relativistic mean-field theory with an effective point coupling between the nucle- 
ons, three-dimensional angular momentum projection is implemented for the first time to project 
out states with designed angular momentum from deformed intrinsic states generated by triaxial 
quadrupole constraints. The same effective parameter set PC-Fl of the effective interaction is used 
for deriving the mean field and the collective Hamiltonian. Pairing correlations are taken into 
account by the BCS method using both monopole forces and zero range (5-forces with strength pa- 
rameters adjusted to experimental even-odd mass differences. The method is applied successfully 
to the isotopes ^"^Mg, ^'^Mg, and ^^Mg. 



PACS numbers: 21.10.-k, 21.10.Re, 21.30.Fe, 21.60.Jz 



I. INTRODUCTION 

Experimental and theoretical studies of nuclei far from 
the /3-stability line are at the forefront of nuclear science. 
Until 1985 [l|], the access to nuclei near the border of 
/3-stability was practically impossible. The advent of ra- 
dioactive ion beams (RIBs) provides a useful tool for 
studying the structure of such unstable nuclei. Hitherto, 
RIBs have already disclosed many structure phenomena 
in exotic nuclei with extreme isospin values, and the next 
generation of radioactive-beam facilities will present new 
exciting opportunities for the study of the nuclear many- 
body systems 0, i, i, 0, i, i, [l3l ■ 

Energy density functional (EDF) theory in nuclear 
physics is nowadays the most important microscopic ap- 
proach for large-scale nuclear structure calculations in 
heavy nuclei and it has been successfully employed for 
the description of nuclei far from /3-stability [ill, fl^ - 
The nuclear EDF is constructed phenomenologically, 
based on the knowledge accumulated within modern self- 
consistent mean- field (SCMF) approaches built upon an 
effective density-dependent two-body interaction. Com- 
pared with the shell model approach [l^, [l3| , EDF func- 
tional are universal in the sense that they can be ap- 
plied to nuclei all over the periodic table. Because of its 
simplicity SCMF approaches have a great advantage in 
particular for the description of heavy exotic nuclei. 

The great success achieved by SCMF theories in the 
description of nuclear properties relies on the fact that 
within these theories the complicated many-body wave 
functions are approximated by a single Slater determi- 
nant. Important many-body correlations are taken into 
account via the mechanism of "spontaneous symmetry 
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breaking" [1^. Examples are the violation of SO (3) ro- 
tational symmetry in deformed nuclei and of U(l) sym- 
metry in gauge space in superfluid nuclei. As a conse- 
quence, such product wave functions are not eigenstates 
of the angular momentum and particle number opera- 
tors. These deficiencies give rise to several serious prob- 
lems in the description of particular nuclear properties, 
as the absence of correlations associated with the symme- 
try restoration, the admixture of low-lying excited states 
into the ground state, difficulties in the connection to 
the laboratory frame for spectroscopic observables, the 
absence of selection rules for transitions, etc. Therefore, 
in order to compare properly with the experimental data, 
one has to go beyond the mean-field approximation. Pro- 
jection methods provide an effective tool to restore the 
spontaneous breaking of symmetries [l^, [13, [H, [l^, HO] ■ 
A suitable linear combination of intrinsic states deformed 
in Euler space or gauge space recover rotational or gauge 
symmetry. Such procedures are known as Angular Mo- 
mentum Projection (AMP) or Particle Number Projec- 
tion (PNP) methods. 

Angular Momentum Projection has been a goal of nu- 
clear physicists for many years. However, due to its nu- 
merical complexity, only in the last ten years it has been 
possible to apply such projection procedure in the con- 
text of SCMF theory with realistic effective forces, for 
example the non- relativistic Skyrme force SLy4 [2l| , the 
Gogny force DIS |2^. [23l or the relativistic point cou- 
pling force PC-Fl 1241251. These investigations have 
shown that the energy gain due to the restoration of ro- 
tational symmetry is of the order of several MeV and it 
has great influence on the topological structure of the nu- 
clear potential energy surface (PES). In these three cases, 
however, axial symmetry in the mean-fields has been im- 
posed from the beginning. Such a restriction simplifies 
the numerical problem considerably, because in this case, 
the integrals over two of the three Euler angles in the 
kernels can be treated analytically and one is left with a 
one-dimensional integration. 
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As illustrated by recent systematic calculations [26l |. 
specific combinations of single-particle orbitals near the 
Fermi surface and the additional binding energy due to 
non-axial degrees of freedom can enhance the tendency 
to form nuclei with triaxial shapes. Several islands of tri- 
axiality have been revealed throughout the nuclear chart. 
The inclusion of triaxiality can dramatically reduce the 
barrier separating prolate and oblate minima, leading to 
structures that are soft or unstable for triaxial distor- 
tions [l^l • Furthermore, the occurrence of triaxiality can 
give rise to many very interesting modes of collective 
motion, which are very different from those of axially 
deformed shape, such as Chiral rotation Wobbling 
motion [2§| and the violation of iiT-selection rules in elec- 
tromagnetic transitions [s^l . 

To describe properly the properties of possible tri- 
axially deformed nuclei and especially to examine the 
role of triaxial deformation in the context of SCMF the- 
ory, it is essential to introduce the 7-degree of free- 
dom at the mean-field level and to perform full three- 
dimensional angular momentum projection (3DAMP). 
In the context of phcnomcnological models with small 
shell model spaces and the corresponding effective inter- 
actions, 3DAMP has already been implemented many 
years ago in Refs. [3l|, [H, IH, [H, [s^. The restora- 
tion of rotational symmetry has been shown to have a 
strong influence on the topological structure of the (/3, 7) 
energy surface for transitional nuclei [33| . In particu- 
lar, the correlations taken into account by 3DAMP arc 
found to have a tendency to lower the potential energy 
in the region of strong triaxial deformations [1^ . In the 
context of energy density functionals, 3DAMP has been 
performed on top of Hartree-Fock (HF) with a simple 
Skyrme-type interaction [stI , or with the full Skyrme en- 
ergy functional [3^ . In both cases cranked wave func- 
tions were projected to approximate a variation after 
projection procedure, but pairing correlations were not 
included. Only very recently, 3DAMP-I-PNP with con- 
figuration mixing has been attempted in the context of 
triaxial Hartree-Fock-Bogoliubov (HFB) theory with the 
full Skyrme energy functional [soj . 

During the past decades, relativistic mean-field (RMF) 
theory, which relies on basic ideas of effective field the- 
ory [iol and of density functional theory [4lj has achieved 
great success in describing many nuclear phenomena for 
both stable and exotic nuclei over the entire nuclear chart 
with a few universal parameters [H, [H, [H, 0, [ill . It 
incorporates many important relativistic effects, such as 
the presence of large Lorcntz scalar and vector fields with 
approximately equal magnitude and opposite sign. This 
leads to a new saturation mechanism via the difference 
between the scalar and vector densities, and naturally to 
the large spin-orbit splitting needed for the understand- 
ing of magic numbers in finite nuclei. Moreover, relativis- 
tic effects are responsible for the efficient description of 
spin observables in medium-energy proton-nucleus scat- 
tering using the relativistic impulse approximation (46j 
and for the existence of approximate pscudospin symme- 



try in nuclear spectra [13, [3 . All these features motivate 
further investigations in the framework of RMF theory 
and new efforts to improve its predictive power. 

The extension of RMF theory for the description of tri- 
axially deformed nuclei was first done decades ago [ioj . 
Later it has been employed in many studies on the effect 
of 7 deformation on nuclear properties [sO, HH, [H, Ull ■ 
7-dcformation plays also an important role in the mean 
field description of rotating nuclei in the framework of 
the cranking model [131: the Coriolis operator violates 
axial symmetry and leads to currents and time-odd com- 
ponents in the intrinsic nuclear fields (55| . All these ap- 
plications of triaxial RMF theory are done on the mean 
field level. A full 3DAMP for such cases is still miss- 
ing and strongly desired, especially for the description of 
transitional nuclei. In this work, we apply for the first 
time 3DAMP to restore rotational symmetry for triaxi- 
ally deformed intrinsic states in the framework of RMF 
theory based on point coupling interactions. 

The paper has been arranged as follows. In Sec. HIl we 
present an outline of the relativistic point coupling model 
that will be used to generate mean-field wave functions 
with triaxial symmetry, and we discuss three-dimensional 
angular momentum projection. The method is applied 
for several isotopes, ^''Mg, ^°Mg, and ^^Mg to check the 
numerical accuracy of the code as well as to present 
several illustrative results in Sec. IIIII Finally, a sum- 
mary and a perspective is given in Sec. IIVI Formulae of 
3DAMP, and details about the calculations of contrac- 
tions and overlaps in the relativistic case are collected in 
the Appendix. 



II. FRAMEWORK 

A. The relativistic mean-field theory with point 
coupling 

A detailed description of RMF theory with point cou- 
pling that will be adopted to generate intrinsic wave func- 
tions can be found in Ref. [56| . In order to present a self- 
contained description of our approach we will give here 
a short outline of the relativistic point coupling model 
used in our applications. 

The elementary building blocks of a RMF theory with 
point coupling vertices are 

i^pom, Oe{i,f}, re {i,7m>75,757m:^a.J: (i) 

where ip is the Dirac spinor field of nucleon, r is the 

isospin vector and P is one of the 4x4 Dirac matri- 
ces. There are ten such building blocks characterized 
by their transformation characteristics in isospin and in 
Minkowski space. We adopt arrows to indicate vectors in 
isospin space and bold types for the space vectors. Greek 
indices /i and v run over the Minkowski indices 0, 1, 2, 
3. 
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A general effective Lagrangian can be written as a 
power series in ipOTip and their derivatives. In present 
work, we start with the following Lagrangian density: 

where the Lagrangian density for free nuclcon reads 

z:''™ - ^(*7m5^ - ^H'- (3) 

The four-fcrmion point coupling term is given by 

1 - - 1 - 

^ = -2"-s(V'V')(^'0) - ■^axsii'f^) ■ (ipfi/j) 

-^"y(V'7M'^)(V'7^V') 

-^aTv{i'fjf,ij) ■ ii^T^'^i'), (4) 

which contains scalar-isoscalar, scalar-isovector, vector- 
isoscalar and vcctor-isovector channels. The medium de- 
pendence of the effective interaction has been taken into 
account by the higher order interaction terms 

, ,(5) 

As in the nonrelativistic Skyrme functional [57[ gradient 

terms are essential. They simulate to some extent the 
effect of finite range of the force: 

-^SrvdMr^^i^) ■ 0'' {^jf^^^^). (6) 

In principle, one could construct many more higher order 
interaction terms, or derivative terms of higher order, but 
in practice only a relatively small set of free parameters 
can be adjusted from the data of ground-state nuclear 
properties. The electromagnetic interaction between pro- 
tons is described as usual 

= -^F'^'-F^.-eVV^^M^, (7) 

where e is the charge unit for protons and it vanishes 

for neutrons. The total Lagrangian density ([2]) contains 
eleven coupling constants as, ay, arv, ckTS, Psi "Js, 
jv, <5s, Sv, Sts and 5tv- The subscripts indicate the 
symmetry of the couplings: S stands for scalar, V for 
vector, and T for isovector, while the symbol refer to the 
additional distinctions: a refers to four-fermion term, S 
to derivative couplings, and P and 7 to the third- and 
fourth-order terms, respectively. 



The pseudoscalar 75 and pseudovector 757^ channels 
do not contribute at the Hartree level due to the parity 
conservation in nuclei and therefore we have neglected 
it in the Lagrangian density From the experience 
of RMF with finite-range (RMF-FR) meson exchange, a 
fit, which includes the isovector-scalar interaction has not 
been found to improve the description of nuclear ground 
state observables. This part of the interaction is therefore 
neglected. Consequently, there are nine free parameters 
in RMF-PC model, which is comparable with those in 
RMF-FR model. 

Using the mean-field approximation and the "no-sea" 
approximation, the operators V'(Cr)i?/' in Eq. ([2]) are re- 
placed by their expectation values and become bilinear 
forms of the Dirac spinor ipk for nucleons 

4>{dru ^ midru]^) = J2 ^iMorUk, (8) 

k 

where i indicates S, V, and TV. The sum ^ runs over 

k 

only positive-energy states with the occupation proba- 
bilities Based on these assumptions, one finds the 
energy density functional for a nuclear system: 

Euf[t,ps,J^,A^] = J d\E{r), (9) 
where the energy density 

£{r) = £'^'"(r) + £-'"'(r) + f ""(r) (10) 
has a kinetic part 

f "^'"(r) = T(r) =Y.vl V4(r) {ap + f3m - m) ^^(r), 

k 

an interaction part 

+ ^J,J^ + '-^ij,fr+'ij,^r (12) 

+^^jTv ■ Utv)^ + —Jtv ■ ^i:iTv)t., 

which contains the local densities and currents 

Ps{r) = ^-y^-0fc(r)V'fe(r), (13a) 

k 

Jt^ir) = Y.'^lMrh^Mr). (13b) 

k 

JTvir) = J2''kMr)T^''Mr)- (13c) 

k 

and an electromagnetic part 
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^24) channel respectively. In the case of the monopolc force 
we have kj, = UfcUfc and 



Minimization of the energy density functional ^ with 
respect to -04, gives rise to the Dirac equation (i.e., Kohn- 
Sham equation) for the single nucleons 



(15) 



The single-particle effective Hamiltonian contains local 
scalar S{r) and vector V^''(r) potentials 

5(r) = Ss, y''(r) =S^+f.S^^, (16) 

where the nucleon scalar-isoscalar E5, vector-isoscalar 

S'' and vector-isovcctor S^y self-energies are given in 
terms of the various densities 



r—n.p 

In the case of a i5-force we use 



fkUkVk 

k>0 



r—n^p 



(20) 



where is the constant pairing strength and the pairing 
tensor K.{r) is given by 



K{r) = -2^ /fcUfcWfe|-0fc(r)|' 



(22) 



fc>0 



asps - 



PsPs ^ 
'lv{fvf 



ISPs 



SsAps, (17a) 
SvAj^ + eA^, (17b) 



(17c) 



The pairing strength parameters Gr in the case of 

monopole pairing and Vr for zero range pairing forces 
are adjusted by fitting the average single-particle pairing 
gap 



For ground state of an even-even nucleus one has time 
reversal symmetry and the space-like components of the 
currents }i in Eq. (|13p and the spatial part of the vector 
potential V(r) in Eq. (fT6|) vanish. Moreover, because of 
charge conservation in nuclei, only the 3rd-component of 
isovector potentials contributes. The Coulomb field 
Aq is determined by Poisson's equation. 

In addition to the self-consistent mean- field potentials, 
for open-shell nuclei, pairing correlations are taken into 
account by the BCS method with a smooth cutoff fac- 
tor fk to simulate the effects of finite-range [S^H^, i.e. 
we have to add to the functional ^ a pairing energy 
depending on the pairing tensor k of the form 



E^^,[k,k*] = J2 fkfk'{k~k\VPP\k'k')K 

kk'>0 

with the smooth cut-off weight factors 



k^k' 



(18) 



fk = 



1 



(19) 



1 + exp[(efc -ep - AEr)/pT] ' 

where ej, is the eigenvalue of the self-consistent single- 
particle field, ei? is the chemical potential determined 
through the constraint on average particle number: 
(<I>|iVr|$) = Nr. The cut-off parameters AE^ and 
Pt ~ AEt/10 are chosen in such a way that = 

fe>o 

Nt + 1.65N^^^, where Nr is the particle number of neu- 
tron or proton. 

In the following calculations we use both a monopole 
force and a density-independent (S-force in the pairing 



(A) 



T,k fkvl 



(23) 



to the experimental odd-even mass difference obtained 

with a five-point formula. 

Moreover, the proper treatment of center of mass 
(cm.) motion has been found very important in the bind- 
ing energy of light nuclei [gO, [U, HI] . We adopt the same 
cm. correction to the total energy after variation, as it 
has been used in adjusting the parameter set PC-Fl [60| . 



2mA ^ 



(24) 



where m is the mass of neutron or proton. A is mass 

number and Pcm = '^f Pi is the total momentum in the 
cm. frame. 

The total energy for the nuclear system becomes 



Etot. = £^df[t, ps,j^, A^] + Ep,i,[K, K*] + EZ"- (25) 

To obtain the potential energy surface (PES), the mass 
quadrupole moment is constrained through the quantities 
920 a-iid q22 , which are related to the triaxial deformation 
parameters [3 and 7 of the Bohr Hamiltonian by 



920 
922 



167r 
"T5" 
32^ 



3 

47r 



Ai?2/3cos7, (26a) 



^ARl^Psin,, 



(26b) 
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where Rq = 1.2A^/^ fm. The total mass quadrupole mo- 
ment q is thus given by 



167r 
5 



q = \l ^\/q^o + '2qh■ 



{27) 



We thus obtain mean field wave functions |$(/3,7)) that 

depend on the deformation parameters f3 and 7. In the 
following we abbreviate the pair of deformation parame- 
ters by a single letter q = (/3, 7). 



B. Three dimensional angular momentum 
projection 

The nuclear mean- field wave function |$) is a prod- 
uct of the solutions of the deformed Dirac equation of 
Eq. (fT5|) and therefore it does not have good angular 
momentum. To obtain the collective energy spectrum 
and wave functions with the good angular momentum 
J, it is crucial to restore the spontaneously broken ro- 
tational symmetry. Especially, for triaxially deformed 
states I $((7)) with the deformation parameters q = 7), 
a full 3D AMP is required. 

The wave function l^'^^g^) in the laboratory frame, that 
is an eigenfunction of and with the eigenvalues 
J( J -I- 1) and M, is obtained by projection [lB| 



10 ^E/""" (9)1^^^^^' 9)' 



(28) 



K 



where a = 1, 2, • • • labels the different collective excited 
states. The basis \JMK,q) functions are not just sim- 
ply Wigner _D-functions as adopted in the classical triax- 
ial rotor model but they are determined microscopically 
from the intrinsic state | by projection using the 
operators P^jj^ 

\JMK,q)=P^jj,Mq)). (29) 
The projector-like operator ^as the form. 



pj _ 

^MK — 



2J + 1 

87r2 



dnDi^ji{n)R{n), 



(30) 



with representing a set of the three Euler angles 
{(j>,9,tp) and the measure dfl = dcf) sin 9 d9 dip. D'l,^j^{Q) 
is the Wigner Z?-function with the rotational operator 
chosen in the notation of Edmonds [6^ as R{^) = 
^vpj^^iejy^iipj^ _ rpj^^ effect of P^jx is extracting from 
the intrinsic state \^{q)) the component with an eigen- 
value K of the angular momentum projection along the 
intrinsic z-axis [isl . [6^ . Since K is not a good quan- 
tum number for a triaxial shape, all these components 
must be mixed, which corresponds to the so-called "K- 
mixing" . Considering the D2 symmetry of triaxial shape 



for even-even nuclei, the sum in Eq. ([28)) is restricted to 
non-negative even values of K. The wave function 1^';;^*^) 
is therefore simplified as [H, 



I* 



JM\ 



K>0 



JMK+,q), 



(31) 



where the angular momentum projected ii'-component, 
\JMK+, q), is given by 

\JMK+,q) = [Pijj, + i-iyP^j_j,mq)). (32) 

The expansion coefficients f^^{q) are determined requir- 
ing that the energy evaluated on 1^^^) stationary with 
respect to fa^*{q)- This condition leads to the general- 
ized eigenvalue equation 

E i'^KK' (9; q) - E-lM^j,, (g; g)}/^^' (q) = 0, (33) 

K'>0 

where the overlap kernels O'j^j^, {q; q) are determined by 



Oi^,{q]q) 



\Oi,j,,{q;q) + {-lf-'Oix_X,{q;q) 



+ i-^yO-k_^,{q;q) + {-iyOi^^,{q;q)], 

(34) 

with d = l,H, and Akk' = 1/[(1 + <5xo)(l + Sk'o)] 
O'KK'iqiq) = / dnDp^,mq)\ORmq)). (35) 



The details about the calculation of overlap functions 
($(q)|O.R(0)|$((j)) win be given in the next section. 

The generalized eigenvalue equation II33D is solved in 
the standard way as discussed in Ref. [l5|. It is accom- 
plished by diagonalizing the norm kernel A/'^^/(g; q) first 



^ M!^l,,{q;q)u'^<'{q) = niui^iq). (36) 
K'>a 

The eigenfunctions u'^ {q) form a complete orthonormal- 
ized set 

Y.u*y^iq)ui^'iq)^5KK', (37a) 

?n 

J2<''iqhif{q)=5r^ra'. (37b) 
K>0 

The non-zero eigenvalues (nj^ 7^ 0) of the matrix 
■^KK' (q-i q) used to build the normalized vectors (i.e. 
the natural states) as 
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1 J 

H^^=J2ui''iq)\JMK+,q), (38) 

V "m _f(->0 

which arc orthogonal and define the "collective" sub- 
space. 

In practice, a cut-ofF x is usually introduced to define 
the non-zero eigenvalues, i.e., n;^ > x- In this work, how- 
ever, we do not need such a cut-off. This is because the 
states with zero eigenvalue in norm matrix have already 
been excluded by constructing the collective wave func- 
tion with the help of D2 symmetry as shown in Ea. (|3ip . 
Of course, if one performs GCM calculations, one cannot 
avoid introducing this cut-off. 

The solution of Eq. ([55]) becomes an usual eigenvalue 
problem, 



(39) 



The reduced matrix element of {Jf , q\\Q2\\Ji, q) is given 
by. 



fJiKi 



KiKf 



tj,M 



(43) 



with j = 2J + 1 and Q2f^ — r^Y2fj,. One can evaluate 
the integration over the Euler angles in the interval [0, tt] 
and multiply with the factor 



1 + i-iy'e-''^^'' + e-'^'" + (-iYe-~^{M+K^)7.'\ ^ (44) 



with the collective Hamiltonian given by the matrix ele- 
ments 

{m\H\m')^ ^ J2 n*^^{q)ni^^,u'S{q). 

V"m'^m' K,K'>0 

(40) 

The solution of Eq. ((39|) determines both the energies 
and the weights fi^{q) of nuclear states l^^f/). 



(41) 



C. Evaluation of electromagnetic transition 
probability 

Once the weights fa^{q) of nuclear collective wave 
function l^'j;^*^) are known, it is straightforward to cal- 
culate all physical observables, such as electromagnetic 
transition probability. Some of them provide a good test 
of the accuracy of symmetry restoration which can be 
used to determine a sufficient number of mesh points in 
the integration over the Euler angles in Eq. ([35|) . More- 
over, through the construction of the collective wave 
function in Eq. (j3ip zero eigenvalues of the norm ker- 
nel have been removed. There are subsequently J/2 + 1 
or (J — l)/2 collective states and rotation energy lev- 
els for the even or odd spin J [35|. These levels will be 
assigned into bands according to their B(E2) transition 
probabilities. 

The B{E2) transition probability from an initial state 
(g, Ji, Ui) to a final state (g, J/, a/) is defined by 



B{E2\ q, J„ai g, J/, a/) 



^j-^\{Jf,q\\Q2U,q){'- 
(42) 



The angular-momentum projection performs a trans- 
formation to the laboratory frame of reference. This 
transformation cannot be inverted and therefore, an in- 
trinsic deformation cannot be unambiguously assigned to 
the projected states. Instead, the comparison between 
theoretical and experimental "deformations" should be 
done directly on the basis of B(E2) values and spectro- 
scopic quadrupole moments Q''^\j,a), 



'I67r / J 2 J , 
e^/^ I (J,g||Q2||J,g)(45) 

\ J J i 



Since the B{E2) values and spectroscopic quadrupole 
moments Q'-'*-' ( J, a) are calculated in full configuration 
space, there is no need to introduce effective charges, 
and hence e denotes the bare value of proton charge. 



D. Evaluation of the overlap integrals 

In the following we evaluate the projected matrix ele- 
ments for general many-body operators O 



2 / I 1 

Oi^, = -r^ / dnD'i;K,mq)\ORmq)) (46) 



8tt' 
2J+1 
87r2 



where, for convenience, we have introduced the following 
notation 



|n).«M, ,47, 



n{VL) 
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with (Ojfi) = 1. The rotational overlap 

n{n) = (0|i?(f7)|0) (48) 

is derived in Eq. (|A28|) of Appendix [X] 

Using the generalized Wick theorem introduced in 
Refs. [6^, m, 1131 the overlap functions (0|O|r2) for arbi- 
trary many-body operators O can be evaluated in terms 
of the mixed densities (|A3p 

Pkim = {0\alak\n), (49a) 
4°(n) = {0\aiak\n), (49b) 
4]{n) ^ {0\alal\n)*. (49c) 

In this way we obtain for instance for a local single par- 
ticle operator Q{r) the projected matrix element 

Q'kk' = J d\Q{r)pij,,{r), (50) 
with the projected density 

P'kK'ir) = J dnDi\,{n)p{r-M)n{a). (51) 

where p(r; fi) is the representation of the mixed density 

(|49ap in r-space given in Eq. (|C7|) 
For the Hamiltonian overlap in Eq. p5|) we find 

H'k,k' = / d\Hij,,{r), (52) 

with 

niK'ir) = J dnDi\,{n)n{r-^)n{n), (53) 

where the mixed energy density has the form 



The Coulomb part of the mixed energy density has the 
form 

^"(-;f^) = £p>;")/^v^^, (56) 

Since the exchange term of Coulomb interaction has 

not been included in the parameterizations of relativis- 
tic mean-field energy density functional, it has been ne- 
glected in the energy kernel as well. 

Because of time reversal invariance the spatial parts of 
the currents jv(r) in Eq. (|13b|) . irvif) in Eq. (|13cp and 
the electromagnetic current jcm('') vanish in the mean 
field calculations. This is no longer true for the mixed 
currents in Eq. (jC10|) . Because of time reversal symme- 
try they are purely imaginary. In the present calculations 
we take into account jy (r; Q) and ]tv{i"i ^) but, for sim- 
plicity, we neglect the small contributions of the gradi- 
ent terms of the mixed spatial currents Ajy(r;ri) and 
AjTy(r;ri) in Eq. and the mixed electromagnetic 

current jom(^; 

The pairing part for the (5-force is given by 

-HT"{r- f^) - - ^ nf*{r- ^W^{r- il), (57) 

where the mixed pairing tensors in coordinate space 

K°^*(r; fi) and K^°(r; fi) are given in Eq. (jC15p . For the 
monopole force we have 

wrm = -GrY. (f^) E '^l'S'(^)' (58) 

fe>0 fc'>0 

where the mixed pairing densities K°i*(f}) and n^^!j.,{^) 

in oscillator space are given in Eq. (|C13[) . 

The cm. correction in Eq. ([24)) is evaluated only 
within the mean field approximation at each value of q. 
The quality of this approximation has not been inves- 
tigated so far. In this case, the contribution from the 
center-of-mass motion to the energy levels of different 
spin is the same at a fixed deformation. 



+n^{r-n)+W^"{r-n). (54) 

The kinetic part 

7i'^'"(r,r!) =T(r;17) (55) 

is given in Eq. (|C11|) . The interaction part W"*(r,17) 

has the same structure as the corresponding energy den- 
sity H™*(r) in Eq. (fT^ . We only have to replace the 
densities p{r) and currents j^{r) by the mixed densi- 
ties p{r]Q) and the mixed currents j'^{r;il) derived in 
Eqs. (jC7[) and (jClOp . This is an ad- hoc procedure that 
is used by analogy to the Hamiltonian case . 



E. Symmetries of the overlap integrals 

The imposed symmetries (I?2 symmetry and time re- 
versal symmetry) in the mean-field calculations give rise 
to symmetries in the overlaps {^{q)\0 R{il)\^{q)) and al- 
low the reduction of the integration intervals for the Eulcr 
angles approximate by a factor of 16 (3^. [35|. 

Specifically, the imposed D2 symmetry reduces the 
integration intervals for the Eulcr angles (0, 9, tjj) in 
Eqs. dSSl) and (gSI) to e [0,7r], 9 € [0,7r], t/. e [0,7r]. 
The symmetries associated with the angles 0, ip for the 
Hamiltonian overlap are summarized as follows: 

{HR{(j),9,il))Y = {HR{^j,9,(j})), (59a) 
{HR{(l>,9,tlj))* = (i7i?(7r-0,6l,7r- V)). (59b) 
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Therefore we have to calculate the Hamiltonian and norm 
overlaps for the Euler angles 0, ip explicitly only in two 
regions: a triangle area with -0 e [O,7r/2],0 e [0, ■0] and 
a square area with ?/' € [7i'/2,7r],0 G [0, 7r/2]. Using the 
above mentioned symmetries we obtain the values in the 
remaining regions. 

For the overlaps of an irreducible tensor operator Ta/j, 
one has the following relationships: 

(fA^i?(^ + 0,0,1/;)) = {-ir{fx^,R{d^,e,^)), (60a) 
(fAMi?(0,^,^ + VO> = {n^R{c^,e,^)), (60b) 
(fA^i?(^-0,0,^-^)) = (-l)^(fA-^iZ(0,0,V)).(6Oc) 

The symmetries associated with 6 are summarized as fol- 
lows: 

(M(0,7r-0,V)) = {HRi^,e,^)), (61a) 
(fA^i?(0,7r- 0,^-7/;)) = {-ir{fx^Ri^,0,^)). (61b) 

Details on the derivation of symmetry properties of 
the overlap integrals can be found in Refs. [lil, and 
in Appendix [P] 

The restoration of broken symmetries in density func- 
tional theory is connected with spurious divergencies, 
which have been observed in connection with number 
projection by the Madrid group in Ref. l69l| and in con- 
nection with the GCM-mcthod in Ref. [70|. Divergen- 
cies have also been noticed in the calculation of overlap 
matrix elements between zero-quasiparticle states and 
two-quasiparticle states in Ref. [7l[. The spurious di- 
vergencies in number projection are connected with level 
crossings and occur in gauge space at the value of the 
gauge angle 99 = 7r/2 for levels with the BCS occupation 
numbers = ^- These poles do not occur in theories 
based on one density independent many-body Hamilto- 
nian, if all the terms in the projected energy are taken 
into account in a consistent way, in particular Fock terms, 
contributions of the Coulomb and spin-orbit potential 
to pairing etc (for details see Ref. [63|). This is obvi- 
ously not the case in most versions of density functional 
theory, as for instance in Skyrme or Gogny functionals 
with fractional density dependence [zl] or for all cases, 
where the effective particle-particle interaction is differ- 
ent from the effective particle-hole interaction. Covariant 
density functional theory, as it is used here, is such a case 
and such poles have been found in connection with num- 
ber pro jection before the variation in relativistic theories 
too [73|. In principle the many-body terms of the point 
coupling Lagrangian in Eq. ([5]) lead to integer powers of 
the density dependence, but the Fock terms are neglected 
and the pairing part of the density functional cannot be 
derived from the same Hamiltonian as the mean field 
part. In fact, most of the successful density functionals 
in the literature have the problem of such poles. They 
cause in particular problems in the case of projection 
before the variation [69l. [t^. In addition, the prescrip- 
tion for the evaluation of mixed energy density in analogy 



with the generalized Wick's theory for Hamiltonian based 
case will also lead finite spurious contributions. 

During the years several recipes have been developed 
to deal with these problems. The most simple method 
to avoid the spurious divergencies is by avoiding the pole 
in the integration over the angles, i.e. by avoiding the 
value If = tt/2 in the case of number projection. Of 
course, this does not help for a very fine integration mesh. 
One therefore has to look for a plateau in the projected 
energy as a function of the number of mesh points. More 
recently a method has been developed in Ref. [1^ where 
the projected energy functional is modified and the terms 
containing the dangerous level crossings and leading to 
finite spurious contributions are removed. 

In the present investigations we have not observed the 
spurious divergencies. In particular we have found con- 
vergence in the number of mesh points (see Figs. [71 [SI 
and [21 of section IIII Bp and therefore the plateau condi- 
tion is fulfilled here. This might be connected to the fact 
that we do not carry out a variation after projection. In 
Ref. [7^ such problems have been observed in the case 
angular momentum projection in systems with cranked 
wave functions and odd particle number. Of course, it 
has to be investigated, whether such divergencies can also 
occur in systems with time reversal invariance. Work in 
this direction is in progress. Moreover, the investigation 
of correction from finite spurious contribution is beyond 
the scope of the present work and will be postponed in 
the future study. 



III. RESULTS AND DISCUSSION 

In this section, we discuss 3DAMP+RMF-PC calcula- 
tions in the nuclei ^'^Mg, ^°Mg and ^^Mg. The intrin- 
sic wave functions that are used in the 3DAMP calcula- 
tion have been obtained as solutions of the self-consistent 
RMF equations constrained on the mass quadrupole mo- 
ments. During minimization, parity, D2 symmetry, and 
time reversal symmetry are imposed. The densities are 
thus symmetric with respect to reflections on the x = 0, 
y = and z = planes. The parameter set chosen for the 
Lagrangian density in Eq. ^ is PC-Fl [s^. The solution 
of the equation of motion psp for the nuclcons is accom- 
plished by an expansion of the Dirac spinors in a set of 
three-dimensional harmonic oscillator basis functions in 
Cartesian coordinates with N^h major shells. The basis 
is chosen to be isotropic, i.e. the oscillator parameters 
are chosen as bx ^ by = ^ bo = \/ h/ mujn in order to 
keep the basis closed under rotations [7y, UM ■ The oscil- 
lator frequency is given by Hujq ~ 4:1A~^/^. The Poisson's 
equation for the electromagnetic field is solved using the 
standard Green function method [53. 
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TABLE L The binding energies Eb (in MeV), ciiarge radii Rc (in fm) calculated by the triaxially (Tri.) deformed and by the 
spherical (Sph.) RMF-PC codes using the parameter set PC-Fl in comparison with the available data. Pairing correlation 
is taken into account by the BCS method with (5-forces. In the triaxial calculations the oscillator shell number is chosen as 
Ng-ti — 12 except for ^'^*Pb with A^^h ~ 14. In spherical RMF calculations, both A^h ~ 12 and A^h ~ 20 are chosen for all nuclei. 
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4.589 
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A. Illustrative examples of mean-field calculations 

To illustrate our triaxial RMF-PC mean-field calcula- 
tion, the total binding energies and charge radii of some 
typical spherical nuclei, adopted for adjusting the PC-Fl 
set, are calculated with triaxially deformed and spherical 
RMF-PC approaches with PC-Fl set. The binding en- 
ergies and charge radii, together with the corresponding 
data available are given in Table [D 

It shows that both the binding energies and the charge 
radii given by the triaxially deformed and spherical RMF- 
PC approaches are in good agreement with the data. The 
tiny differences in the binding energies by these two ap- 
proaches are due to the different numerical algorithm. 
Here, we have to point out that the binding energies of 
^"^Ca and ^^Ni with N — Z are relatively poorly repro- 
duced with a difference of about 2-3 MeV which cannot 
be cured simply by increasing the shell number A^sh and it 
may be ascribed to the missing of proton-neutron pairing 
correlations in the present calculations. 



B. Convergence check of three-dimensional angular 
momentum projection 

In Fig.[l] we show the mean- field binding energy curves 
for ^'^Mg as functions of the mass quadrupole moment q 
(<Z22 = 0) defined in Eq. ([27|) . calculated by the triaxial 
RMF-PC approach with the parameter set PC-Fl. The 
four different energy curves correspond to the calcula- 
tions with A^sh = 6, 8, 10, and 12 major oscillator shells 
respectively. It shows that A"sh = 8 is sufficient to obtain 
a reasonably converged mean-field binding energy curve 
for ^^Mg. Pairing correlations have been taken into ac- 
count by the BCS method with monopolc pairing forces. 
The pairing strength parameters Gr are determined sep- 
arately for neutrons and protons by adjusting the pair- 
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FIG. 1: (Color online) Binding energy curves for ^*Mg, cal- 
culated by the constrained self-consistent triaxial relativistic 
mean-field approach in a three-dimensional harmonic oscilla- 
tor basis with major shells Ash ~ 6, 8, 10, and 12 respectively. 



ing gaps of the mean-field ground state to the odd-even 
mass difference as obtained with a five-point formula. 
The pairing strength parameters G„ = 34.6/^ MeV and 
Gp = 33.75/^ MeV determined in this way have been 
kept fixed throughout the constraint calculations. 

In Fig. [5] we plot the pairing gaps of neutrons and pro- 
tons in ■^^Mg as functions of the quadrupole moment q 
(922 ~ 0) together with the corresponding energy curve. 
The total energy shows a prolate deformed minimum in 
the energy curve at g = 1.04 with Etot. = -193.57 MeV. 
This figure indicates clearly that the pairing gap changes 
considerably with the deformation refiecting the changes 
in the single particle level density. Obviously the min- 
imum in the energy corresponds to a rather low level 
density p^ . 

For an axially symmetric intrinsic state, the norm over- 
lap in Eq. (|M2)) can be calculated analytically using the 
Gaussian Overlap Approximation (GOA) [66l. pTSj: 
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FIG. 2: (Color online) Upper panel: Pairing gaps A^^„/p 
of neutron (dash line) and proton (solid line). Lower panel; 
energy curve for ^''Mg in a triaxial RMF-PC+BCS calculation 
with a constant pairing strength Gt , determined by fitting the 
ground state gaps At to the odd-even mass difference. 



n(g,g;O,0,O)«exp[--(j2)sin2 0], 



(62) 



which turns out to be an excellent approximation and 

thus provides a very useful test of the numerical proce- 
dure used in angular momentum projection [23 |. 

Fig. [3] displays the norm overlaps n{q, q; 0, 9, 0) as fimc- 
tions of the Euler angle 9 for several different axially 
deformed intrinsic states of ^^Mg. It shows that the 
3D AMP calculated values of the function n{q; 9) are in 
good agreement with those given by the GOA approxi- 
mation. 

For triaxially well-deformed intrinsic states, the norm 
overlap has been derived approximately in Refs. [tsI. Fzoj: 



« cxp[--(j2)02 _^ (cos(0 -f V') " 

(63) 



+ -{Jx)9{sm.4)- si-a-i))] 



In our calculation for ^^Mg, the third term in the ex- 
ponential vanishes because of time reversal invariance 
{Jx) = 0. The norm overlaps n{q,q;4)/Tp) are given in 
Fig. [3] as functions of the Euler angles (j) and tp for sev- 
eral triaxially deformed intrinsic states of ^^Mg. It is 
found that the norm overlaps n{q,q;(f>/ip) oscillate in an 
exact 3D AMP calculation as functions of </> and ip with a 
period of T = 180°. The approximate formula Eq. ([M]) 
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FIG. 3: (Color online) Left panel: A comparison between the 
norm overlaps n{q, q; 9) as functions of the Euler angle 6 for 
several intrinsic states of ^^Mg obtained by a 3D AMP calcu- 
lation (open circle) and the GOA formula (solid curve) . Right 
panel: A comparison between the norm overlaps n{q,q\(f)/'ip) 
as functions of the Euler angles cj) and i/' for several different 
intrinsic states of ^*Mg obtained by a 3DAMP calculation 
(open circle) and the GOA formula (solid curve), where the 
Eulers angle are = 0.02 and (f) = 0.02 radian with 9 = 0.04 
radian. 
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FIG. 4: (Color online) The mixed densities p{r;q,q;(f>,9,^p) 
in x-z plane with /3 = 0.55,7 = 0°. The Euler angle 9 has 
the value of 1.23°, 61.28°, 118.72°, and 178.77° respectively 
keeping <f) = tp — 0. 



is obviously valid only in the interval 0° to 90°. In or- 
der to obtain the approximate results in the interval be- 
tween 90° and 180° we use symmetry around the angle 
= 90°. In this case. Fig. [3] shows that the Gaussian 
overlap approximation can roughly reproduce the results 
obtained by the exact 3DAMP calculations. Moreover, 
as expected, the larger the 7 deformation of the intrin- 
sic state is, the larger is the amplitude of the oscillating 
norm overlaps n{q] 4>/ip)- 

To describe the collective motion of nuclei in the con- 
text of energy density functional theory, one should de- 
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FIG. 5: (Color online) The different terms of the Hamiltonian 
overlap for (3 = 0.55 and 7 = 0° as functions of the Euler angle 
6 normalized to S = 0°. 
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FIG. 6: (Color online) Upper panel: Hamiltonian overlap with 
/3 = 0.48, 7 = 0° as a function of the Euler angle 6. Lower 
panel: Hamiltonian overlap with /3 = 0.40, 7 = 26° as a func- 
tion of the Euler angles cj) and %p. The values with 6 = Q or 
(j>{ or xp) — are chosen as zero. 



tcrminc the corresponding collective Hamiltonian. In the 
3DAMP+RMF-PC approach, the matrix elements of col- 
lective Hamiltonian are constructed in Eq. ([34]) in terms 
of the Hamiltonian overlaps, which have their standard 
functional form but depend upon the mixed densities and 
currents. 

In Fig. m we plot the mixed nucleon densities 
p(r; q, q; 0, 0, jp) in the x-z plane derived from the mean- 
field state with /? = 0.55,7 = 0° for </) = V = 0° 
and for various Euler angles 9 ^ 1.23°, 61.28°, 118.72°, 
and 178.77°. It is obvious that the reflection symme- 
tries with respect to the planes x = 0, y = and 
z = present in the mean-field densities are violated in 
the corresponding mixed densities. Moreover, we show 
in Fig. [5] the various terms in the Hamiltonian over- 
lap h{qa,qa',^) = a{0\H\n)a resulting from the four- 
fermion coupling term, the current contributions, the 



Coulomb term, the derivative term, the kinetic term, the 
higher order term and the pairing term as functions of 
the Euler angle 9 for the mean-field state at the point 
f3 = 0.55,7 = 0°- The energy surface is normalized to 
6* = 0, i.e. Ah{qa, qa] ^) = Klaila] ^) - h{qa, qa, 0). We 
find that the current contributions and the Coulomb term 
in the Hamiltonian overlap change mildly with the rota- 
tion angle 9 and thus they have only small contributions 
to the collective Hamiltonian. On the contrary, the four- 
fermion coupling term, the pairing term and the higher 
order term are sensitive to the Euler angle 9 and play a 
dominant role in the collective Hamiltonian. 

In Fig. [SI we display the total Hamiltonian overlap for 
the axially deformed mean-field state with (3 = 0.48, 7 = 
0° and the triaxially deformed mean-field state with (3 = 
0.40,7 = 26° as functions of the Euler angle 0, or the 
Euler angles (j) and ip. It shows that both for the axially 
deformed shape and the triaxially deformed shape, the 
Hamiltonian overlaps, behaving like the norm overlaps, 
oscillate with the period T — 180° in the Euler angle 6*, 
4>, or 

A N-point Gaussian-Legendre quadrature is used for 
integration over the Euler angles 0, 9 and ip in the cal- 
culations of the norm kernel M^j^, and the Hamiltonian 
kernel Ti-xx'- The calculation of the Hamiltonian over- 
lap at each mesh point of the Euler angles is very time 
consuming. Therefore, besides the utilization of symme- 
tries in overlaps, it is essential to make a careful check 
of the convergence for the number of mesh points. The 
projected energy and the B{E2) transition probability 
are good observables for this purpose. 

In Fig. [71 we plot the projected energy of first O"*" state 
obtained from the mean-field states with (3 = 0.55, 7 = 0° 
and (3 = 0.21,7 = 0° for ^"^Mg, and the corresponding 
B{E2 [: 2+ — > 0+) transition probabilities as fimctions 
of the number of mesh points ne for the Euler angle 9. 
The projected energy of the 0"*" state from the mean-field 
states with (3 = 0.55,7 = 28° and the B{E2 [: 2+ -> 
0^) transition probability, as functions of the number of 
mesh points (or n^) for the Euler angle (p (or ijj) are 
shown in Fig. [51 where 9, (j> and ip have values between 
and TT. We find that in order to achieve a precision 
of 0.001% for Eo+ and 0.1% for B{E2 : 2+ ^ 0+) the 
total number of mesh points for the Euler angles in the 
intervals cp e [0,7:], 9 e [0,t:], tp € [0,7r] should fulfil the 
relation: x Ng x > 6 x 6 x 6. 

According to the uncertainty principle A J • AQ ~ h, 
we need a large number for meshpoints in the Euler 
angels for higher values of the spin. In Fig. [51 we 
show the projected energies of 2+,4+ and 6+ obtained 
from the mean-field states with (3 ~ 0.47,7 = 17° and 
(3 = 0.52,7 = 0° as functions of the number of mesh 
points 71,0, n.0 and ng. We find that it is possible with 
X Ne X > 12 X 14 X 12, to achieve a precision of 
0.001% in the energy of a projected state with angular 
momentum up to J = 6 in the ground state band. In the 
following calculations we use such large numbers of mesh 
points in the Euler angles. 
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FIG. 7: (Color online) The projected energy of the lowest 0"*" state derived from mean-field states with /3 = 0.55, 7 = and 
f3 = 0.21,7 — 0° for the nucleus '^^Mg, and the B{E2 [: 2+ — > 0"'") transition probability, as functions of the number of mesh 
points ne for the Euler angle 9. 
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FIG. 8: (Color online) The projected energy of the first 0''' 
state obtained from the mean-field state with /3 = 0.55, 7 = 
28° for ^*Mg, and the B{E2 j: 2+ 0+) transition probabil- 
ity, as functions of the number of mesh points (or n^) for 
the Euler angles (p (or -0). 



FIG. 9: (Color online) The projected energies of 2+,4+ and 
6^ states from the mean-field states with /3 = 0.47, 7 = 17° 
and {3 = 0.52, 7 = 0° for ^*Mg as functions of the number of 
mesh points (or n^) or ng. Eq is the converged energy of 
a state with spin J. 



Since the states with very small occupation probabili- 
ties give negligible contributions to the kernels, as usual, 
we introduce a cut-off parameter C, which divides the 
Dirac space into an occupied part and an unoccupied 
part (see Eq. (jA14[) . The states with vl < C will be ex- 
cluded in the calculation of the overlaps. In Fig. (TU] we 
show the projected energy of the first 0+ state and the 
B(E2 .[■. 2+ — > O"*") transition probability projected from 
the mean-field state with /3 = 0.55, 7 = 0° as functions of 
the cut-off parameter It shows that ( should be chosen 
as C < 10^^ m order to get a precision of 0.01% for Eq+ 
and of 0.00001% for the B{E2 i: 2+ ^ 0+) value. Using 
the cut-off C reduces the computational effort ( about 
80% of total computer time for Ngh = 8) in the calcu- 
lations of the norm overlap and the matrix elements of 
mixed densities and pairing tensors considerably, espe- 



cially for the cases of large iVsh, small particle number 
and weak pairing, where most single particle levels of the 
Dirac basis have nearly zero occupation probabilities. 



C. Tests of three-dimensional angular momentum 
projection 

1. Application to an axially deformed shape 

To illustrate the validity of our newly-developed 
3DAMP-|-RMF-PC-hBCS code, we first apply it to the 
axially deformed case, where a IDAMP calculation is 
possible. The projected ^ = 0+, 2+, 4+, 6+, and 
8"'' potential energy curves of '^^Mg have already been 
calculated with IDAMP+RMF-PC-f BCS approach in 



13 



LU 
CO 




FIG. 10: (Color online) The projected energy of 0^ state from 
the mean- field solution with /3 = 0.55, 7 = 0° for ^*Mg, and 
the B{E2 [: 2+ 0"*") transition probability, as functions of 
cut-off C, in Dirac space ). 



Ref. [2J]. To make a comparison, we perform the 
same calculations within the 3DAMP-f RMF-PC+BCS 
approach. The numerical techniques are the same as 
those of Ref. [23|. We find that our newly-developed 
3DAMP4-RMF-PC-hBCS code can reproduce the results 
given by IDAMP-hRMF-PC-hBCS approach. 

Furthermore, following Ref. [1^, we first test the 
3DAMP+RMF-PC approach for an axially deformed 
shape, which allows two distinct orientations in the in- 
trinsic frame: the symmetry axis can cither parallel to 
the 2r-axis or it can be perpendicular to it. 
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FIG. 11: (Color online) The excitation spectra and B{E2) val- 
ues projected from the axially deformed mean-field states with 
(3 = 0.44, 7 = 180° and (3 = 0.44, 7 = 60° respectively. The 
first and second columns show the unique band with K — 
and the unique band with Tf-mixing. The last four columns 
show the decomposition into A'-components when the sym- 
metry axis is chosen perpendicular to the z axis, i.e. the 
K — 0,2,4,6 bands respectively. 

In Fig. [TI] we show the excitation spectra and B{E2) 
values for ^''Mg projected from the axially deformed 
mean-field states with f3 = 0.44,7 = 180° and (3 — 
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FIG. 12: (Color online) Upper panel: The probabilities W'' 
of finding the component with given spin J. Lower panel: the 
probabilities of finding the component with given spin J 
and the projection K along the z-axis in the mean-field state 
with the deformation parameters /3 = 0.44, 7 — 60°. 



0.44,7 = 60° respectively. For 7 = 180°, z-axis is along 
the symmetry axis, and therefore only one pure K — 
band can be found. All other if-components have zero 
norm. While for 7 = 60°, one can show that the pure 
K = Q state is transformed into a multiplet of states 
with K ranging between and J. Such phenomena 
can be seen more clearly from the probabilities W'^ = 
^^($((7)|Pj^^|(f>((7)) of finding a component with given 

spin J and the probabilities = {^{q)\PKK\^{<l)) /W'' 
of finding a component with given spin J as well as given 
projection K along 2-axis. These probabilities are shown 
in Fig. [HI 

In principle, the transformed wave functions differ only 
by an unobservable phase and the energies of projected 
states as well as the electromagnetic transition probabil- 
ities should be identical. This provides us an excellent 
test of the numerical accuracy of the projection scheme 
in the code. Fig. [11] shows that for the low spin states, 
e.g., O'*", 2+, the projected energies and B{E2) values are 
exactly the same. As angular momentum increases, the 
difference increases to a largest value (~ 0.4%) in the 
B{E2 j: 6+ 4+), which could be reduced with more 
mesh points in the Euler angles. 



2. Application to a triaxially deformed shapes 

The excitation energies and B{E2) values for ^^Mg 
projected from the triaxially deformed mean-field states 
with (3 = 0.61,7 = 10.3°; P = 0.61,7 = 109.7° and 
P = 0.61, 7 = 130.3° are presented in Fig.[Tl AU the ex- 
citation energies are arranged into bands according to the 
B{E2) values. These three intrinsic states correspond to 
the same nuclear shape with three different orientations 
in the intrinsic frame. The projected energy and the 
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FIG. 13: (Color online) The excitation spectra and B{E2) values projected from triaxially deformed mean-field states with 
P = 0.61, 7 = 10.3°; (3 = 0.61, 7 = 109.7°; and p = 0.61, 7 = 130.3° respectively. 



electromagnetic transition probability do not depend on 
the orientation of the nucleus and therefore, in principle, 
the predicted values should be the same as illustrated in 
Fig. [m It shows that the projected energies and B{E2) 
values in these cases are in good agreement with each 
other. However, small differences in the B{E2) values 
appear and increase with angular momentum. Except 
for the B{E2 : 6+ 4+) in the K = 2 band, the differ- 
ence is smaller than 1%. This indicates that more mesh 
points in the Euler angles are necessary to provide a bet- 
ter description of the B{E2 : 6+ ^ 
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3. Dispersion of particle numbers 
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FIG. 14: (Color online) The probabilities W'' in the mean- 
field states with f3 = 0.61,7 = 10.3°, = 0.61,7 = 109.7° 
and f3 — 0.61, 7 = 130.3° respectively. 



The decomposition of a triaxial mean-field state into 
components with different J-values in the laboratory 
frame should also be independent on its orientation in 
the intrinsic frame. In Fig. 1141 we show almost the same 
probabilities W"' of different spin states in these cases. 
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FIG. 15: (Color online) The average neutron (filled circle) and 
proton (open circle) numbers of angular momentum projected 
states with J < 6 from axially deformed intrinsic states of 
2*Mg with P = 0.2,0.3,0.4. 

Although the mean-field intrinsic states are obtained 
with the constraint on the right average particle num- 
ber, it cannot guarantee the right particle number in the 
angular momentum projected states. In order to make 
up this flaw, in principle, one has to perform PNP cal- 
culation. The study with both PNP and 3D AMP in the 
context of GCM has only been attempted based on a 
Skyrme EDF theory [3^. Such kind of study based on a 
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covariant EDF theory is still extremely time-consuming. 
As the first step, in this work, neither the exact pro- 
jection on particle numbers N and Z, nor a constrain on 
the average number of particle in the angular-momentum 
projected states is performed. Therefore, it is essential 
to know the dispersion of particle numbers within a rota- 
tional band. In Fig. 1151 we plot the average neutron and 
proton numbers of angular momentum projected states 
with J < 6 from axially deformed intrinsic states of ^^Mg 
with j3 = 0.2, 0.3, 0.4. It shows that the error in average 
particle number of projected states with J < 6 is within 
0.5%. 



D. Examples of three-dimensional angular 
momentum projection 

1. Application to '^'^ Mg 



TABLE II: Reduced E2 transition probabilities from states 
Jf to states J] in ^''Mg. The minimum of the projected 
J = 2 PES is used for the calculation of intrinsic wave 
function. The experimental data for the excitation ener- 
gies iSx [in units of MeV] and E2 transition probabilities [in 
units of e'^fm'*] are taken from most recent available sources, 
le^fm'' = 'iTT{lf{1.2A^^'-^)-^ W.u. = 0.243 W.u. for ^^Mg. 
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The 3DAMP+RMF-PC approach has been used in 
Ref. [ill to describe the PES in the /3-7 plane for the low- 
est = 0"*" state and for the first excited = 2+ in the 
nucleus ^^Mg. There is no pronounced minimum with an 
obvious 7-dcformation in the PES for the 0"*" state, which 
is in disagreement with the results of Ref. [3^, where 



a 3DAMP+PNP calculation based on a non-relativistic 
Skyrme HFB functional shows a pronounced triaxial min- 
imum with 13 ~ 0.6 and 7 = 16°. Keeping in mind that 
we found strong pairing gaps in our mean-field calcula- 
tions, an additional number projection is not expected to 
change this result. A possible reason for this difference is 
the fact that different energy functionals are used in these 
two calculations. A minimum with /3 w 0.55, 7 ~ 10° has 
been found on the PES of the first excited 2+ state. To 
construct the excitation spectrum and to calculate B{E2) 
transitions, one should in principle perform a GCM con- 
figuration mixing calculation on top of three-dimensional 
angular momentum projection, or choose the minimum 
of the different J projected PES as basis. However, such 
kind of calculations are beyond our present study. In- 
stead, wc use the minimum of the projected J = 2 PES 
as basis to calculate the experimentally observed excited 
energy levels and the B{E2) transition probabilities in 
^'^Mg using Eqs. ^ and This is the only way to 

obtain K = 2 bands in our calculation. The details about 
the B{E2) transition probabilities in ^''Mg are given in 
Tab. nil It shows that the predicted intraband B{E2) 
values arc systematically smaller than the data, while 
the interband B{E2) values are systematically overesti- 
mated. It indicates that the amplitude of "K-mixing" is 
too strong in our calculations. 
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FIG. 16: (Color online) The probabilities W'' and the prob- 
abilities in mean-field states with (3 = 0.55 as functions 
of triaxial deformation 7. 



In order to understand the effect of 7-deformation on 
the amplitude of angular momentum mixing, it is useful 
to investigate the individual components forming the in- 
trinsic state |$((7)), i.e. the J- and K-xmKvag. In Fig. [TBI 
we present the probabilities W"^ and the probabilities 
in the mean-field states with (3 — 0.55 as functions of 
the triaxiality parameter 7, ranging between 0° and 60°. 
It is noted that each 7-deformation in this range corre- 
sponds to a definite shape uniquely. Fig. [11] shows that 
the J-mixing remains practically constant with changes 
in the 7-deformation, while the amount of i^-mixing in- 
creases considerably with increasing triaxiality. This in- 
dicates that the underestimated intraband B{E2) values 
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and the overestimated interband B{E2) in the low- lying 
excited states of ^^Mg as shown in Tab. [II] are due to the 
large 7-deformation in the intrinsic state. 

To illustrate the effect of 7-deformation on the B{E2 [) 
values, we plot in Fig. [T7] the intraband B{E2 J,) transi- 
tion probabilities for 2^ — > 0^, 4]*" ^ 2^ and 6^ i'l in 
the ground state band projected from mean-field states 
with P = 0.55 as functions of the 7-dcformation. Ob- 
viously the intraband B{E2 |) values increase when 7 
approaches 0° or 60°. It indicates that a configuration 
mixing calculation (GCM) within a generator coordinate 
method might be very important to understand the ob- 
served B{E2) values. Alternatively, calculating B{E2) 
value using the minima of each J projected PES might 
also improve the results. Moreover, we note here that 
in contrast to the B(E2 I) values for the 4^ 2f and 
61 4j^ transitions with obvious minima at 7 ~ 10°, 




B{E2 i) values for 2+ ^ Oj^ 
with 7, ranging from 72 e^fm'' 



FIG. 18: (Color online) The potential energy surfaces of 
mean-field state and projected 0^ states in the /3-7 plane ob- 
tained by triaxial RMF-PC-I-BCS calculations for the nucleus 
^°Mg. The contour lines are separated by 0.5 MeV. 



obtained at MSU and at GANIL using the method of 
intermediate-energy Coulomb excitation are 295(26) e^ 



changes only moderately 

to 88 e^fm'* which is con- ^™ and 435(58) e fm* [8J|, respectively. However, 



sistent with the data 86.4 ± 1.6 e^fm^ of Rcf. 
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FIG. 17: (Color online) B{E2 [) transition probabilities for 
O;;'', 4+ — > 2i and &'l 4+ in the ground state band, 
projected from the mean-field state with jS = 0.55 in the nu- 
cleus ^*Mg, as functions of the triaxial deformation 7. 



2. Application to ^° Mg 

The evolution of shell structure and appearance of new 
magic numbers in neutron-rich nuclei has become one of 
the main topics in recent investigations of nuclear struc- 
ture physics. Especially, the erosion of the neutron magic 
numbers iV = 20 and 28 and the occurrence of well- 
deformed prolate deformed structures in such magic or 
close-to-magic nuclei are presently in the focus of several 
investigations. 

There is much controversy about the deformation of 
the ground state in the nucleus ^°Mg. Experimen- 
tally, this deformation is determined by measuring the 
B{E2; 0+ 2+) transition probability. The values 



the most recent measurement performed at CERN re- 
sults in 241(31) e^ fm^ [8^1, which is lower than those 
extracted in previous measurements performed at inter- 
mediate energies. Therefore it is very interesting to study 
this problem theoretically within the present approach. 

In Fig. [18] we plot the potential energy surfaces of 
mean-field states and projected 0"*" states in the P-j plane 
for the nucleus •^°Mg. The intrinsic states are calculated 
in the triaxial RMF-PC-I-BCS approach using monopole 
pairing forces with G„ = 24.4/A, Gp = 29.7/A, adjusted 
to the experimental odd-even mass differences. We find 
that the mean-field potential energy surface is very soft 
against f3 in the spherical region. It is hard to recognize a 
minimum. The energy surface projected on the 0+-state 
has, however, a pronounced axially symmetric minimum. 

Fig. [19] shows axially symmetric results for the nu- 
cleus states in ^°Mg. The corresponding potential en- 
ergy curves of the intrinsic states and of the projected 
J'^ = 0^',2+,4+,6+ states are given as functions of the 
quadrupole moment q ((722 = 0). The intrinsic deformed 
states are obtained in the RMF-PC-I-BCS approach us- 
ing either a monopole pairing forces or a zero range S- 
type pairing forces. We find that the projected curves 
for the 0"*" state have in both cases an obvious minimum 
at /3 ~ 0.25. The energy differences between the min- 
imum and the spherical shape are 3.87 MeV (BCS-G) 
and 3.69 McV (BCS-(5) respectively. The correspond- 
ing B{E2 : 0+ 2+) values are 194.8 e^fm'' and 194.6 
e^fm'*, respectively. Both of them are somewhat smaller 
than the data. 



3. Application to '^'^ Mg 

For the nucleus '^^Mg, a much lower excitation energy 
of 0.885 MeV was measured for the first 2"'"-state [8g 
and a large deformation with /? ~ 0.51 has been inferred 
from the measured B(E2: 0+ ^ 2+) value (454 ± 78 



e fm ) [87[. Therefore this nucleus has drawn much at- 
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FIG. 19: (Color online) Potential energy curves of the pro- 
jected J" = 0+, 2+, 4+, 6+ states in ^°Mg, as functions of the 
quadrupole moment q. The intrinsic deformed states are ob- 
tained by RMF-PC-f BCS calculations with both monopole 
pairing forces (left panel) and (5-type pairing forces (right 
panel) . The pairing strength parameters Vt for the zero range 
pairing forces are adjusted the experimental pairing gap as 
discussed in Eq. (f23)) . 



FIG. 20: (Color online) The neutron (left panel) and proton 
(right panel) single-particle levels for '^'^Mg, as functions of 
the quadrupole deformation j3. The levels with positive (neg- 
ative) parity are shown with solid (dashed) lines. The levels 
belonging to the f-r/2 orbit are plotted with red dashed lines. 
The fermi energies for neutrons and protons are plotted with 
blue dotted lines. 



tention in studies with self-consistent approaches. Cor- 
rections from the angular momentum projection and con- 
figuration mixing are found to be essential to reproduce 
the large deformed ground state of '^^Mgin the HFB 
approach with the Gogny force DIS [H, [8§|. However, 
similar non-relativistic calculation with Skyrme-type the 
Sly4 force [9^ and relativistic calculation with the PC-Fl 
force fail to reproduce the data [2^ . Therefore, it is inter- 
esting to revisit this problem within our 3DAMP-f RMF- 
PC approach. 

Fig. [20] displays the neutron and proton RMF- 
PC-I-BCS single-particle energy levels for '^^Mg as func- 
tions of the quadrupole deformation /3. The pairing 
strength parameters are Gn = 26.78/A and Gp = 
32.25/v4 for the monopole pairing force. They are ob- 
tained by adjusting the gaps at the spherical minimum 
(the ground state of the mean-field calculation) to the 
experimental odd-even mass difference with a five-point 
formula. In the self-consistent calculations we find a col- 
lapse of proton pairing for the range 0.45 < (3 < 0.75. 

The potential energy curves of the projected J'^ ~ 
0+,2+,4+,6+ states in ^'^Mg are plotted in Fig. [21] as 
functions of the quadrupole moment q ((722 = 0). The in- 
trinsic deformed states are obtained from RMF-PC+BCS 
calculations with monopole forces and 5- forces. The pair- 
ing strengthes Vr are adjusted to the odd-even mass dif- 
ference. 

At the mean-field level, a shoulder of only 1.8 MeV 
above the spherical minimum has been found in the 
present calculations with pairing strength parameters ad- 
justed to odd-even mass differences. This value is close to 
the prediction of 1.9 McV for the shoulder by the HFB 




FIG. 21: (Color online) Potential energy curves of the pro- 
jected J" = 0+, 2+, 4+, 6+ states in ^^Mg, as functions of the 
quadrupole moment q. The intrinsic deformed states are ob- 
tained by RMF-PC-I-BCS calculations with monopole forces 
(left panel) and (5-forces (right panel). The pairing strength 
parameters Vt for the zero range pairing forces are adjusted 
the experimental pairing gap as discussed in Eq. (|23|l . 



approach with the Gogny force [23| . but much smaller 
than the value of 3.5 MeV predicted by the RMF-PC 
model with ^pairing forces taken from the parameter set 
PC-Fl set [ill. In Fig. [H] we show various RMF cal- 
culations for this shoulder with the parameter sets PC- 
Fl [ill, PKl [m and NL3 [H. Pairing correlations are 
taken into account by the BCS method with a monopole 
pairing force (BCS-G) or a (5- force (BCS-i5). In all cases 
the pairing strength parameters are adjusted to the odd- 
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FIG. 22: (Color online) The energy curves as functions of 
mass quadruple moment g (7 = 0°) for "^^Mg, calculated with 
different RMF parameterizations. 



60 60 




FIG. 23: (Color online) The potential energy surfaces of 
mean-field theory (left panel) and of angular momentum pro- 
jection J — after the variation (right panel) in the /3-7 plane 
obtained by triaxial RMF-PC+BCS calculations for ^^Mg. 
The quadrupole deformation of the minimum in potential en- 
ergy surface of 0^ state is /3 ~ 0.6, 7 ~ 10°. The contour lines 
are separated by 0.5 MeV. 



even mass difference except the case labeled by "BCS-i5*" 
where Vt has bee taken from the PC-Fl set |24| . 

We find that the energy curves in RMF calculations do 
not depend on too much on the effective interactions but 
rather strongly on the strength of the pairing force. All 
the calculations with a pairing strength adjusted to the 
experimental pairing gaps give a lower shoulder, while the 
calculation with a iS-pairing forces taken from the PC-Fl 
set produce a higher shoulder with a stiffer energy surface 
against quadrupole deformation (3. Similar phenomena 
have also been found in Skyrme-Hartree-Fock-|-BCS cal- 
culations [92|. As a consequence, one will obtain different 
predictions for the deformation of ground state for dif- 
ferent pairing correlations. More detailed investigations 
concerning this question are in progress. 

In Fig. [23] we examine the potential energy surface in 
the /?-7 plane for ^^Mg. Triaxial RMF-PC-fBCS cal- 
culations with a monopole pairing force (left panel) are 
compared with angular momentum projection on J = 0. 
We observe that considering the 7-degree of freedom 
one can expect considerably enlarged ground state de- 
formations. The quadrupole deformation of the mini- 
mum in the angular projected 0"'' PES is found to be at 
(3 ~ 0.6,7 — 10°i based on which, the predicted energy 



of the 2+ state is E = 1.21 MeV and the predicted B(E2: 
0+ 2+) value is 573.5 e^fm''. It has to be pointed 
out that the PES of 0+ state is very /3-soft in the re- 
gion 0.3 < (3 < 0.7. Based on the intrinsic state with 
quadrupole deformations (3 = 0.3,7 = 0, the AMP pre- 
dicted an energy of 2+ is £■ 3.39 MeV and a B(E2: 
0+ _> 2+) value of 250.2 e^fm^. This indicates clearly 
that the generator coordinate method based on 3DAMP 
approach becomes necessary for a full understanding of 
the properties of ^^Mg. 



IV. SUMMARY AND PERSPECTIVE 

In this paper, a full three-dimensional angular mo- 
mentum projection on top of a triaxial relativistic mean- 
field calculation has been implemented for the first time. 
The underlying Lagrangian is a point coupling model 
and pairing correlations are taken into account both 
with a monopole force and a (5-force. Convergence has 
been checked and the validity of this newly-developed 
approach has been illustrated by applying it to the de- 
scription of the low-lying excited states in several Mg 
isotopes. 

For ^^Mg no pronounced minimum with obvious triax- 
ial deformation has been found on the potential energy 
surface of the O"*" state. A minimum with /3 « 0.55, 7 « 
10° has been found on the PES of the first 2+ state. Using 
this minimum as a basis for the projection the experimen- 
tally observed excitation energies and B(E2) transition 
probabilities can be qualitatively reproduced. However, 
the predicted spacing between the levels is overestimated 
in this approach. 

For '^"Mg, the projected energy surface of the 0"^ state 
has a obvious minimum with (3 ~ 0.25. The energy differ- 
ences between the minimum and the spherical shape are 
3.87 MeV (BCS-G) and 3.69 MeV {BCS-S) respectively. 
The corresponding B{E2 : O"*" 2+) are respectively 
194.8 e^fm" and 194.6 e^fm*. 

For '^^Mg, we note that the calculations with ad- 
justed pairing strength parameters produce always a 
lower shoulder in the mean-field energy curve, which is, 
together with the triaxial degree of freedom, essential to 
reproduce the large deformed ground state. Moreover, 
the mean-field and the projected 0+ potential energy sur- 
faces of ■^^Mg have been found to be very 7-soft in the 
region of small deformations and /3-soft in the neighbor- 
hood of its minimum. 

These investigations indicate that, besides triaxiality, 
the effects of pairing correlations and shape fluctuations 
should be treated more carefully in the description of 
low-lying excited states of exotic nuclei. Work in this 
direction is in progress. 

Finally, we would like to point out that the pairing 
strength parameters of protons and neutrons in PC-Fl 
are adjusted to the pairing gaps of the nuclei: ^^^Xe, 
^'^"'Sm, ^^^Sn, ^^°Sn and ^^^Sn respectively. However, 
pairing strength parameter obtained in this way might 
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not be well-justified in the region of light nuclei. We have 
found in the present study for ^"Mg: fp = 1.04, /„ = 1.19 
and for ^'^Mg: fp = 1.51, /„ = 1.09, where fr is the 
ratio of the pairing strength parameters of the adjusted 
delta-pairing and the standard PC-Fl delta-pairing. It 
indicates that a better parameterizations of the energy 
density functional is required for the description of light 
nuclei. 



tions {^{qa)\OR{n)\^{qb)) for an arbitrary many-body 
operator O can be expressed in terms of the mixed den- 
sities and mixed pairing tensors 



Pkiiqa,qb;^) 
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APPENDIX A: EVALUATION OF 
CONTRACTIONS AND OVERLAPS 

The contractions and overlaps have been derived in de- 
tail in Ref. [2l|, where, however, the rotation matrix is 
assumed to be real from the beginning. This is no longer 
the case for a three-dimensional angular momentum pro- 
jection. In addition these earlier investigations were done 
only for nonrelativistic density functionals. Therefore, 
we derive here in a similar way the general formulae 
of the contractions and overlaps suitable for the three- 
dimensional relativistic case and used in the present nu- 
merical applications. 



1. Determination of the generalized contractions 

In the following we derive formulae of generalized con- 
tractions {^{qa)\OR{n)\^{qb)} connecting different in- 
trinsic states. Such formulae can be applied directly in 
future Generator Coordinate (GCM) calculations with 
3DAMP as well. 

For convenience, we introduce the following notation 



In order to derive expressions for these mixed densities we 
consider the fact that the quasiparticle operators (a, a^) 
and f^/?^) are connected by a Bogoliubov transforma- 
tion [13| 



a \ ( Ut VM / /3 



(A4) 



On the other hand, the quasiparticle operators a, are 

related to the particle operators a,a^ by a Bogoliubov 
transformation. 



(A5) 



In a similar way the quasiparticle operators (3, fH'' are re- 
lated to the particle operators b, by 





(A6) 



Assuming that the operators a, and b, 6^ are related 
by a rotation as [73 



R{n) 
R*{n) 




(A7) 



one finds for the particle operators a, and the quasi- 
particle operators /?, the relation 



^mqa)\, \n)b 



mmqb)) 



(Al) 



mqa)\Rin)mqb)) 

The quasiparticle vacua |0)a and \Q)b arc defined by the 

corresponding quasiparticle operators ak and (3k respec- 
tively, 

a{0\ai=0, Pk\^)b^O. (A2) 
According to the generalized Wick theorem the contrac- 



(A8) 



with the coefficients t/fc(ri), Vf,(rJ) given by 

Ubin) = Rin)Ub, Vbifi) = R*in)Vb. (A9) 

Combining Eq. (|A5P and Eq. (jASj we obtain the matrices 
U and V in Eq. (g4| 
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= UlR{n)Ub + V^R*in)Vb, (AlOa) 

yt = ulRin)vc + v^R*{n)u^, (Aiob) 

that relates the quasiparticle operators /?, /Jt and a, 
and the quasiparticle vacua a(0| and |il)f, in Eq. (jAip . 
With the help of generalized Wick's theorem [63|, one 
finds the contraction 

a{0\al3^\n)b = V-\ (All) 

and in combination with Eqs. (|A2[) . (|A5p and (|A8p . the 

elements of the mixed density and the mixed pairing ten- 
sors of Eq. (|A3|) are obtained as 

PM{qa,qb;n) = [K(r!)[U^]-VJ]H, (Al2a) 
4°{qa.,qt;n) = [VCmiV^r'U^U, (A12b) 

4liqa,qt;n) = [u^mv^rwn^k- (A12c) 



2. Restriction to the occupied space 

In practical three-dimensional applications the matri- 
ces U, V, V etc. have the very large dimension of the 
oscillator basis. In fact most of the high-lying eigen- 
states of the Dirac equation are not occupied and there- 
fore they do not contribute to the overlap integrals. In 
order to reduce the computational effort it is therefore 
of great importance to eliminate these high-lying eigen- 
states in the Dirac basis, where the mean field wave func- 
tion has the form of a BCS wave function. The procedure 
discussed in the following is, however, not restricted to 
RMF-fBCS calculations used in this investigation. In 



general Hartree-Bogoliubov theory one can apply similar 
formulae in the canonical basis p^ . [9^ where an arbi- 
trary Hartree-Bogoliubov wave function has BCS form. 
In this basis the intrinsic states are characterized 

by the special Bogoliubov-Valatin transformation of the 
form 

uJ"" »V 'A. ,A13) 

\0 Uk) \-Vk j 

Here Uk,Vk are real positive numbers and the phase has 
been chosen as uj. = Uk, v^. = —Vk, where k is the time 
reversed state of k. Since unoccupied states with ~ 
have no contribution to overlap and contractions, one can 
eliminate these states to simplify the calculation [2l|, It^, 
[q^ I . As usual, one can introduce a cut-off C to divide the 
full Dirac space into two parts: an occupied part with 
> and an unoccupied part with < C and the 
matrices U and V in Eqs. (|A5p and (jA6p have the form 

V=\ \^U=\ , i?= . (A14) 

\o oy ly \roi Roo) 

The matrix R is related to the occupied states only. In 
this case the matrix IJ-'" in Eq. (|A10a|) becomes 

^ / U^R*U; + VJRV* UlRl,\ ^^^^^ 
and its inverse has the form 



D- 



rJq {r^ 



D-W^{Rn-'Rl 

RUR^-'u^D-^im 



-1 dT 
^01 



(A16) 



where the matrix D is defined as, 

D = UjiR'y'U^ + V;fRVb*. (A17) 
In the general case of GCM calculations where qa ^ qb 

the BCS-space of the wave function \^{qa)) is different 
from the BCS-space of the wave function |$(qf,)) and 
therefore the cut-off procedure can lead to occupied sub- 
spaces and to matrices Va and Vb with different dimen- 
sions and rectangular matrices R and Z), which cannot 
be inverted. In such cases appropriate cut-off parameters 
Ca and have to be chosen, such that the matrix D stays 



a square matrix. 

The elements of the mixed density in Eq. (jA12[) are 

/ RWD-^Vj' \ 
Pkl{qa,qb;n)^ - - • (A18) 

The matrices i?oi and i?io connect the occupied space 

with the unoccupied space by rotation. We neglect these 
matrices in the mixed densities and pairing tensors, be- 
cause they are usually very small, i.e. we restrict our- 
selves to the occupied space in the further calculations. 
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In this space we obtain the elements of the mixed density 
and the mixed pairing tensors as: 

PM{qa,qb;^) = [RV:D-W;[]ku (Al9a) 

4°((7a,gb;^^) = [RVCD-WJ]u, (A19b) 
Rfi{qa,q,;n) = [i?*[7*D-ivf];,. (A19c) 

This shows that we finaUy have to invert only the matrix 
D in the occupied subspace. The explicit expressions for 
the matrix elements of D in Eq. (jA17P are 

Du = ul{lF)llu\ + vlRlivl (A20a) 
D^l = ul{R')-lu\ + vlR*^vl (A20b) 

where the indices fc, I run over the states with non- 
vanishing occupation numbers. Using the time reversal 



properties of the rotational operator, one finds the fol- 
lowing relations: 

Du = ~DIj, Du^Dli. (A21) 

3. Determination of tlie overlaps 

The norm overlap has already been derived in Ref . [67l | , 

=±VdetU. (A22) 
After some calculations we obtain 



= uTR*ut -f vTrv^ = 



1 v^Rvcu;-^Rl, 



U'^R^-^U; + V^RVC 




1 ■ 



(A23) 



and using det R = 1 we find that the norm overlap in 
Eq. (jA22p is simply a product of two determinants of 
much smaller dimension: 

det U = det 75 det R. (A24) 

The phase of the overlap in Eq. (jA22p remains open. 

Neergard and Wiist pointed out that the phase prob- 
lem of the norm overlap could be avoided by rewriting the 
norm overlap into the following form 

k>Q 

= iU^'kOll^'^ + oi), (A25) 

k>0 l>0 

where Uk and u'f, are the Bogoliubov-Valatin transforma- 
tions coefficients in (|A13[) for the intrinsic states |(f>(ga)) 
and |$((7b)) respectively. The product Y[i>o runs over the 
pairwise degenerate eigenvalues c; of the matrix M [t^ 

M{qa,qb;n)^ z,{n)zl with z = v*u*-^ (A26) 

In the canonical basis the matrix Z is reduced to 2 x 2- 
matrices of the form 




(A27) 



where k runs only over states with v"^ > In cases, 

where some of the numbers Uk vanish one can, in analogy 
to Eq. (jAl4p . reduce the space intro three subspaces of 
fully occupied state (w^ = 1), partially occupied states 
{0 < vl < 1) and empty states {v^ = 0). Finally, the 
norm overlap can be evaluated according to Eq. (|A25[) 
by diagonalizing the matrix M. This method is certainly 
rather complicated. It turns out that we do not need to 
apply it in the present applications based on time reversal 
symmetric wave functions \^{q)). The norm-overlap is 1 
for = 0, it stays real and positive for all values of the 
Euler angles il = {</>, 0, -0) and therefore we have for the 
norm overlap 



mqa)\RmMqb)) = VdetDdeti?. (A28) 



APPENDIX B: REPRESENTATION OF 
ROTATIONS IN THE DIRAC BASIS 

In our calculations, the single-particle wave functions 
ipk are Dirac spinors. For the solution of the Dirac equa- 
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tion the large and small components /(r, s) and g{r,s) 
of a Dirac spinor are expanded in terms of the eigen- 
functions of a three-dimensional harmonic oscillator in 
Cartesian coordinates 14911 



iY.9nk\n) 



(Bl) 



where Xtkif) is the isospin part. The harmonic oscillator 
basis states \n) = \nx,ny,nz,ns) with simplex Us ~ +i 
and the time reversed states |n) with simplex Ug ^ —i 
are defined by 



\n) = (l>nAx)'l^ny{y)(l>nAz) 



i^y I 1 



(B2a) 



In) = 0„Jx)</.„^(2/)0„Jz)i-^ r |,(B2b) 



where the phase factor i"" is consistent with the triaxial 
self-consistent symmetries and leads to real matrix ele- 
ments for Dirac equation [i^ [o^, [o^l • 

The matrix elements of the rotation operator R{Q) in 
Eq. (|A14p in the Dirac basis are derived from the rep- 
resentation of this operator in the harmonic oscillator 
basis (|B2| by 



d^r^plir,qa)RmMr,qb) 



= Y.fnkiqa)fn'i{qb){n\Rm\n') (B3) 

n.n' 

+ J2 9^k{qa)gn'iiqb){n\Ri^)\n'). 



The rotation matrices {ni\R{Cl)\n2) in the cartesian basis 
have been derived using the method of generating func- 
tions in Ref. [H]. In present work, however, we adopt 
a simple method to evaluate these matrix elements by 
transforming from the cartesian basis to the spherical 
oscillator basis given by |m) = \nrljm) with 



m 



(B5) 



where the expansion coefficients Fkm and Gkm can be 
obtained with the help of relation in Eq. (|B4|) . 



F, 



rnk 



Gfri 



771 k 



k{m\n) 



(B6) 



The matrix elements of R in Eq. (jB3p are subsequently 
given by 

Rkiiqa,qb;n) = J2 P^^kiqa)F,n'iiqb){ni\Rmm') 

mm' 

+ E G*rr^k{Qa)Gr^'l{qb){m\Rmrh'), 
fhifi' 

(B7) 

where the matrix 

{m\R{n)\m') = S^^n'M^n'DlU^) (B8) 

is diagonal in the quantum numbers Ur, I, j and is simply 

given by the Wigner D-function. We use Condon-Shortly 
notation for the spherical harmonics Yi^ {9, 4>) [i^l . With 
the time reversal operator 

\m) = f\nrl3ni) = (-l)'+J'-">^,;j - m), (B9) 

one finds the expansion coefficients of the Dirac spinor 
for the time reversed state, 

(BIO) 

where |— m) = \nrlj—m) and where k is the time reversed 

state of k. With these relations, the matrix element, Rf.j 
can be easily calculated. 

Moreover, according to the time reversal properties of 
the rotational operator i?(ri), one immediately finds: 



Ru 



kl 



-R 



kV 



Rki = R 



kl- 



(Bll) 



mi ms 

(B4) 

where C;'™ , ,„ is the Clebsch-Gordon coefficient. The 

transformation coefficients {nrlmi\nxnynz) are given in 
Therefore, the large and small compo- 
nents /(r, s) and g(r, s) of ipk can be rewritten in terms 
of the eigenfunctions of spherical harmonic oscillator as. 



APPENDIX C: MIXED DENSITIES IN 
COORDINATE SPACE 

In a point coupling model with a local interaction of 
zero range the overlap integrals for the Hamiltonian are 
most easily evaluated in coordinate space. We therefore 
have to calculate the mixed local densities and currents in 
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r-space. Expressing the Dirac spinors in terms of spher- 
ical harmonic oscillator states 



iGk{r,a) 

with the large and small components 



Ffc(r,cr) = ^^;„fc$„(r,cr), 

rn 



(CI) 

(C2a) 
(C2b) 



and the spherical oscillator functions 

mirris 

Here Xm ^^e spin part. 

According to Eqs. (|A13|) and (|A19a[) we obtain the rel- 
ativistic mixed single-particle density matrix in the har- 
monic oscillator basis, 



rmm' 



-+ ^ 
mm' 



Pfhrh 



G^{n)V^*D-^V'''^G'''^ 



(C4a) 
(C4b) 
(C4c) 
(C4d) 



where the rotated large and small components of Dirac 
spinor, F„ik and G^k are given by 

F.mk{n) = Y,R^"n'{n)F^,k, (C5a) 

Grnkm = ^i?™™-(r!)G™-fc. (C5b) 

fh' 

For an arbitrary one-body operator O, such as the mul- 
tipole moment operator Ta^, the corresponding overlap 
is determined by the mixed density, 

rain' 

m-jn' P^n'iji 

{qa, qb;n). {C6) 

rhfh' 

Finally we obtain for the mixed densities in coordinate 
space 

p{r;qa,qf:n) = ^ P++ (r) |$,„(r)) (C7) 

mm' 

±IZPmm'(*m'(^)l*™(r)), 



where the lower sign holds for the scalar density ps in 
Eq. (|13a|) and the upper sign for the vector density pv 
in Eq. (|13bp . The rotation operator i?(f2) does not com- 
mute with the reflections on the x = 0, y = 0, and z = 
planes. Therefore one has to extend the coordinate rep- 
resentation of the mixed density p{r;qa,qb;^) from 1/8 
to 1/2 of the full space, leaving only parity and isospin 
projection as good quantum numbers. 

Considering the fact that the time reversal operation 
T commutes with spatial rotations R(fl) and time re- 
versal invariance of the quasiparticle vacua: T\0)a = 
\Q)a,T\^l)b = \^)b; one finds that the contributions from 
spin up and down to the mixed density p(r; qa, qb', ^) are 
complex conjugate to each other, 

p{r,(T;qa,qb;^) = p*{r,-a;qa,qb;^), (C8) 
where the relation 



-2aa 



(C9) 



has been used. This shows that the mixed densities 

p(r; qa, qb', ^) in coordinate space, summed over the spin 
index a are real. 

Moreover, there are non-vanishing mixed currents 
j(r; (7a, (7f,; f2) with matrix elements of the same form as 
the densities. 

j{r;qa,qb;n) = -i pt^^, {^n^' {r)\cT\^M) 

mm' 

+^Y.P^''m'{'^m'ir)\rT\^Ur))- (CIO) 

mm' 

Since the total wave functions |$(t", q)) are invariant un- 
der time reversal, these real part of these currents van- 
ishes. 

The mixed kinetic energy in Eq. (j55p is given by 
T{r-qa,qb;n) = - X]P™n'(*m'(r)|cr. V|$,„(r)) 

mm' 

+ T.Prnm'i'^^n'ir)\a-V\^rr.{r)) (Cll) 

mm' 

- m[pv{r; qa, qb; ^) - Ps{r; qa, qb; ^)]- 
Using time reversal invariance and 

a{0\aj,ak\n)b = -a{0\aka-k\n)l, (C12a) 
a{0\alal\n)l = -a{Q\alal\n)l, (C12b) 

we obtain for the mixed pairing tensor in Dirac-space 



'^kk 



^01 



kk 



"kk 



(C13a) 
(C13b) 
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and we find in analogy to Eq. (|C4p for the mixed pairing 
tensors in oscillator space 



10++ ^ 
mm' 

10— ^ 
mm' 

,^01++ _ 



01 — 



~lr7aT 





mm 


(C14a) 




mm' 


(C14b) 


fJapaT 


- * 

- mm 


(C14c) 




- * 

- mm' 


(C14d) 



The many-body Hamiltonian H is rotational invariant, 

which leads to together with orthogonality to the follow- 
ing symmetry relations for the Hamiltonian overlap 

{HR{<j>,9,^)} = (D4) 
{HR{d,,e,^)Y = {HR{~^,-e,~^)), (D5) 
{HR{<j>,9,^)y = {HR{4',0,^)). (D6) 



With the help of relation: e'^'^'-'^e '^•^ae "''^ 
gets 



and in coordinate space 



K ir]qa,qb;^) 



E <?,t+<fnv(r,a)a>„,(r,a) (C15a) 

im' .a 

E 4?m'"$m,(r-,a)ci>,-„(r,a), (C15b) 

im' .a 

E ^m^'^nAr,cTWmir,<j) (C15c) 

E '^"L^KAr,c7)^Ur,a) (C15d) 



In this investigation, GCM and configuration mixing 
is not taken into account. Therefore we have |0)a = |0)h 
and only diagonal contractions with qa — qt ~ q- 



APPENDIX D: SYMMETRIES IN OVERLAPS 

1. Symmetries associated with (p and tp 

The Z?2 symmetry and time reversal symmetry have 
been imposed in the mean-field calculation, which leads 
to the mean-field state |<I>((7)) invariant under the follow- 
ing transformations. 



{HR{(l>,e,ij)y = {HRi-i^, -9, -<!>)) 

= (M(7r -(/), 61, 7r-V')>, (D7) 
which can also be derived from the reality condition: 

{HR{(j>,d,ij)r = {HRi-<j>,e,~^j)) 

= (i7i?(7r-0,6l,7r- V)). (D8) 

In a similar way we can derive symmetries of the over- 
laps with {Tx^R{Q)) . Since Tx/i is not rotational invari- 
ant, the overlaps with the Euler angles (/>,')/' in regions 
[0,7r] and [tt, 27r] are related by the following relations, 

(f A^7?(^ + 0, ^, V^)) - (-l)^(fA^i?(0,0,V^)),(D9a) 
{fxM,e,Tr + ^P)) = {fx^R{<j),0,^)). (D9b) 



The tensor T^^ is transformed under e 



-ITT J-j. 



(DIO) 



= k = x,y,z. 



(Dl) 



which gives rise to the symmetry: 
{fx^Ri(t>,9,i^)) = {-l)\fx-^R{7T-(l),e,7r-i;)). (Dll) 



It reduces the integration intervals for the Euler angles 
(</),6',i/') in Eqs. ([#1) and gSl) to S [0,7r], G [0,7r], 
ijj G [0, tt]. The Hamiltonian kernel Hj^j^, and the norm 
kernel Nj^j^, are simplified as 

Stt^ Jo Jo Jo 
x{dR{<i},9,^))Di\,{<j},9,^), (D2) 

where O = 1,H and the factor Tkk' = 1 + e~''^^ + 
g-iK TV j^^-i{K+K )7r Pm-thcrmore, the rotation operator 
R{4', 9, Tp) is transformed as 

e-*''-^=i?(0, 9, ij)e'''^^ = -9, -tp). (D3) 



2. Symmetries associated with 9 

Since the mean- field state is invariant under the 

transformation e"'^", 

{HR{(I), n-9,i')) = (iJ e"^'^= er^'^-'y e'^'^y e*'^'^^ ) 
= {HR{,p,-9,^^)) 
= {HR{cP,0,n-4,))*. (D12) 

On the other hand, the group elements in the group D2 
obey the relation: e"-^" = e^'^'^^e*'^'^^. 
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(HR{<j>,TT~0,i^)) = {He 



icl>J^ iirjj, iirj^ p-iSJy itfiJ^ 



(D13) 



This shows that the Hamihonian overlap is reaL With 
the help of the relati 
finds the symmetry, 



the help of the relation: e''^'^=TA^e~*^'^*- = (-l)^f\^, one 



These symmetries of the hamiltonian overlap integrals 

simplify the calculations considerably by reducing the 
necessary interval, where the overlap integrals have to 
be calculated from [0, tt] to [0, 7r/2]. 
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